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We study conformally-flat initial data for an arbitrary number of spinning black holes with ex- 
act analytic solutions to the momentum constraints constructed from a linear combination of the 
classical Bowen-York and conformal Kerr extrinsic curvatures. The solution leading to the largest 
intrinsic spin, relative to the ADM mass of the spacetime eg = S/M^uMy is a superposition with 
relative weights of A = 0.783 for conformal Kerr and (1 — A) = 0.217 for Bowen-York. In addition, 
we measure the spin relative to the initial horizon mass Mhq , and find that the quantity x = S/Mfj^ 
reaches a maximum of ■^^^'^ — 0.9856 for A = 0.753. After equilibration, the final black-hole spin 
should lie in the interval 0.9324 < Xflnai < 0.9856. We perform full numerical evolutions to compute 
the energy radiated and the final horizon mass and spin. We find that the black hole settles to a 
final spin of Xfinai = 0.935 when A = 0.783. We also study the evolution of the apparent horizon 
structure of this maximal black hole in detail. 
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I. INTRODUCTION 

With the breakthroughs in the numerical techniques to 
evolve black-hole binaries (BHBs) [iHSl : Numerical Rela- 
tivity (NR) has become a very important tool to explore 
highly-dynamical and nonlinear predictions of General 
Relativity. NR can novif be used to evolve generic black- 
hole binaries and accurately compute gravitational radia- 
tion from such systems. However, these computations are 
very costly when one explores various corners of the BHB 
parameter space, such as extreme mass ratios and nearly 
maximally spinning black holes. Yet those cases are of 
great astrophysical interest [1]. Recently, BHBs with a 
mass ratio of 1:100 were successfully simulated [SHZ]j and 
equal-mass BHBs with spin parameters up to at least 
0.97 of the maximum spin were evolved for 25 orbits ^. 
These latter simulations used the first order generalized 
harmonic formalism and excision of the black-hole (BH) 
interiors in a pseudospectral evolution scheme. 

Simulations of highly-spinning black holes are impor- 
tant for understanding astrophysical black holes. Impor- 
tant spin effects have been found in the merger of BHBs, 
such as the hangup effect |^, large recoil velocities of 
up to 5000km/s [TU], and their statistical consequences 
for galaxy merger |llj . Recent observational evidence for 
highly recoiling black holes was reviewed in Ref. |12) . 

Many of the groups in numerical relativity have 
adopted the so called, moving punctures formalism [31[3], 
where the interiors of the black holes are evolved without 
imposing internal boundary conditions. Since this ap- 
proach usually starts from conformally flat puncture |13j 
initial conditions, and since there is a proof (under the as- 
sumption of axisymmctry and that such a slicing should 
vary smoothly with the spin parameter a) of the nonex- 
istence of conformally flat slices of a Kerr black hole [M] , 
it is generally believed that these data cannot represent 



highly- spinning black holes. However, the above state- 
ment is not quantitative, nor does it take into account 
the other "half" of the initial data, i.e. the extrinsic 
curvature. Even if we could find slices that make the 
Kerr metric conformally flat, they would presumably not 
be maximal, i.e. X 7^ 0. This means that the Hamil- 
tonian and momentum constraints would couple, mak- 
ing the latter nonlinear and superposition of solutions no 
longer valid. 

The first and most widely used initial extrinsic cur- 
vature data for spinning black holes is the Bowen-York 
(BY) solution [15^ to the momentunr constraints. This 
assumes a conformally flat 3-metric and a longitudinal, 
trace-free extrinsic curvature. Cook and York [1^] stud- 
ied intrinsic spins of these data [normalized to the square 
of Arnowitt-Deser-Misner (ADM) mass [T7] (ADM) mass 
of the spacetime], eg = S'/A'/^Qjyj , and found that eg ap- 
proaches a maximum value of e'J'g'^ — 0.928. Further in- 
creases of the spin for these data leads to a corresponding 
increase in the ADM mass such that the intrinsic spin is 
unchanged. Later, in Ref. [TB], it was shown that the con- 
formal extrinsic curvature of a Kerr BH was also an exact 
solution to the conformally flat momentum constraints. 
We denote such initial data by "cKerr" to distinguish 
them from the true Kerr initial data. Those solutions 
have a maximum intrinsic spin of ^^'^^^„ = 0.932. This 
raises the possibility that, by choosing a different extrin- 
sic curvature, one can reach still higher maximum values 
of the intrinsic spin, while keeping the conformal flatness 
condition for the initial 3-inetric. In this paper we will 
explore this possibility explicitly. Assuming the optimal 
extrinsic curvature is still close to that of a Kerr BH, 
we will study a parametrization of deviations from the 
conformal Kerr extrinsic curvature proportional to the 
difference between the conformal Kerr and BY extrinsic 
curvatures. This plays the role of a sourceless extrinsic 



curvature wave. We also perform high resolution evolu- 
tions that allow us to shed some light on the transition 
from the near maximally spinning conformally flat initial 
data to the final submaximal spinning Kerr black hole. 
In particular we track the evolution of trapped surfaces 
and apparent horizons. 



II. INITIAL VALUE PROBLEM FOR BLACK 
HOLES 

We assume a conformally flat 3-metric for the initial 
configuration 

Ir] = V'^'t^J , (1) 

where 7^- is a conformal 3-metric and the 3-Ricci tensor 



Ri 



0. The physical extrinsic curvature is given by 



K,.j = V~^^y, and K'^ = ip-'^^K'-K The indices are 
raised and lowered with the conformal 3-metric 7^^. V; 
is the covariant derivative with respect to the conformal 
3-metric 7^^. 

Assuming the maximal slicing K = Ki' = Q, the mo- 
mentum constraint becomes, 



V^K.j = STTJi 



(2) 



where ji is a conformal matter momentum density, and 
the physical matter momentum density is given by ji = 

The Hamiltonian constraint the leads to 



Ail) + \k,jk'''ilj-^ ^0, 



(3) 



where, for conformally flat initial data, A is the ordinary 
Laplacian operator. Taking i/; = f + mp/{2r) + u(x), 
where nip is a parameter, r = \x\, and u{x) is continuous 
function (in fact, C^) leads to unique solutions of the 
Hamiltonian constraint. 

To solve the momentum constraint in Eq. (pi) , we split 
the conformal extrinsic curvature Kij into a transverse- 
traceless and longitudinal part 19J of the form 



K,, = (Ly),j+T, 



I] , 



where the operator L is given by 



{l^V\ 



2V(,l^,)--,5,,V,y\ 



and Tij is a symmetric, transverse-traceless tensor, 

V^ T,j = , T,'=0. 
Here Vi must be a solution of 



(4) 
(5) 

(6) 

(7) 



The spinning BY extrinsic curvature ^15., is obtained 
from the vector field 



V 



1 



Ajk 



UjSk , 



(8) 



where n^ = x^/r in the Cartesian coordinates {x, y, z}, 
5' is a spin vector, and the Levi-Civita symbol is defined 
as e^y^ — e^.,,, = I. When we include r = in this 



analysis, the momentum constraint becomes 



V'K. 



BY 



ATre,jkS''V^S'^^\x'^ 



(9) 



If we do not add an additional Tij part, the spinning BY 
solution is a purely longitudinal solution. 

To solve for these initial data numerically, we 
use the TwoPunctures thorn [5D] and the Cac- 
tus/Carpet/EinsteinToolkit infrastructure pTti24] . 
We restrict our analysis here to single BH spacetimes. We 
solve the Hamiltonian constraint for i/j using the TwoP- 
unctures thorn by setting the mass and spin of the 
second puncture to zero. 

We use AHFinderDirect |22 to locate apparent 
horizons ( AHs) . We measure the magnitude of the hori- 
zon spin using the Isolated Horizon algorithm detailed 
in Ref. j26J. Note that once we have the horizon spin, 
we can calculate the horizon mass via the Christodoulou 
formula 



rriH 



ml. + Sj,/{4mfJ, 



(10) 



where m[„ = •\/A/(167r) and A is the surface area of the 
horizon, and Sh is the spin angular momentum of the 
BH (in units of A'P). We denote the horizon mass on 
the initial slice and the final equilibrated horizon mass 
by Mhq and Mh^ , respectively. 



III. GENERALIZATION OF 
CONFORMALLY-FLAT INITIAL-DATA 



In Ref. [18] it was shown that 
2 



K,. 



u 



(11) 



is a solution to the momentum constraints, where cos 9 — 
S'^rii and S*' ~ S^/S, provided that the coordinates are 
(quasi-) isotropic (i.e. the conformal metric is flat, or 
conformally Kerr in quasi-isotropic coordinates). The 
(otherwise arbitrary) scalar function lu gives the angu- 
lar momentum of the BH [15] . That is, the spin angular 
momentum Ji is given by 



J. = 



1 
16^ 

_^ 
4 



e,,kf {x^K''"'-x''K^"')d^S.r, 



(ioie ^ t:) - Lu{e ^ 0)\ 



(12) 



assuming -0 — > 1 in the limit ?' — > 00. 

The spinning BY initial data can be recovered by tak- 
ing 



(jJby = —S (cos'^ 9 — 3 cos ( 



(13) 



On the other hand, Eq. (Ill is hnear in to, and we can where 



therefore take an arbitrary superposition J2i Ci Wi, where 
Ci are constants and each uii yields a solution to the vac- 
uum momentum constraint. In particular, we are inter- 
ested in 



/3 = 



16r^S 

3 



= 16 r'' -HMk'' + a"* + 24 r^AfK^- 8 r^a^ 
2 MK^a^ 



5uj 



,BY 



Mk a^ sin* cos 



(14) 



+8 tMk - 2 AfK-'a^ - 8 r A/ko 
-f 32 r^Afx + 16 aV^ cos^ 6* . 



(18) 



where a = S'wavc/^K, S = (rsL^ + a^ cos^ 0) and tbl = 
r[l -h (Mk -f a)/(2r)] [1 + (Mk - a)/(2r)]. Here we denote 
the spin parameter by Swavc rather than S because it is 
in a transverse-traceless "wave-like" part of the extrinsic 
curvature (see below), and does not contribute to the 
ADM angular momentum (unlike S). Therefore, it need 
not match the BH's spin. 



As a way to study the effects of changing the extrinsic 
curvature on the maximum intrinsic spins, we set the 
conformal extrinsic curvature to 



K,,. 






ASK, 



tj , 



(19) 



and vary the parameters A, Mk, and the spin parameter 
S'wavo = AfKa used to construct SKij — {Kfj — ^ij^)- ^^ 
Replacing ui in Eq. (11) by Suj yields a solution to all cases the spin used to construct Kf^ remains fixed 
the momentum constraints with zero spin, which we will at 5* = Af^. 
denote by SKij = Kfj — K^^ . We note that SKij is 
transverse-traceless, i.e., \7^SKij = everywhere, be- 
cause it is constructed from the difference of two ex- 
trinsic curvatures each leading to the same source term 
ji in Eq. (|2|). Finally, the extrinsic curvature given by 
K,, = Kf^ 



ASKij solves the momentum constraints, 
with non-zero spin, for any fixed value of A. Taking A = 1 
yields a pure conformal Kerr extrinsic curvature (pro- 
vided that one uses the same spin parameter 5'wave = S 
to construct K^^ and SKij), while A = corresponds 

to pure BY curvature. Again, we note that SKij does 
not contribute to the spin, so in principle, one can use a 
different spin parameter to construct K^^ and SKij. 



IV. ANALYSIS OF SINGLE SPINNING BLACK 
HOLE INITIAL DATA 



The extrinsic curvature in Eq. (11) can be written as 



K- 



- [{ae - {S'"'n„-,far)eM{i S^ n^ n^) 



-fa^(5™n„)efe£(,S''=n^4)] 
where the coefficients as and a^ are given by 



BY 



as; 



a; 



BY 



0, 



(15) 



(16) 



for the BY case, and for Suj we have 



Sae ^ -16MKa3[(-40rMKa^ + 5a4 + 5AfK'^ 

-10 Mk^o^ + 80 r'^ - 24 r^a^ + 40 rMy^ 
+160r3A/K + 120r2AfK^) cos^ 
- (2r + A/k + af (2r + A/k - af 
+48 a^r'^ cos'' 0\ / (3^ , 

Sar = 32 Mk a^ (2 r -f Mk -t- a) (2 r + A/k - a) 

x{a^-MK^+Ar^) /p^ , (17) 



A. ADM mass dependence 

To estimate the final (equilibrated) the intrinsic spin 
for a large number of configurations, we measure the 
ADM mass of a single spinning BH spacetime with spin 
parameter S = M^ and puncture mass rup = 0.02M. We 
choose this puncture mass so that it has a negligible effect 
on the total ADM mass, while still leading to a solvable 
initial data problem (taking rUp — would lead to sin- 
gularities at the puncture). For the sequences presented 
here, we only perform numerical evolutions to calculate 
the equilibrated spin on a select set of configurations. For 
the rest, we note that equilibrated value of the intrinsic 
spin Sh /Mfj (as measured using the equilibrated hori- 
zon mass) must be larger than Sh/MJ^y^m^ which holds 
true because these axisymmetric data can radiate mass 
but not angular momentum. Since the total mass ra- 
diated is smaU {SM/M < 0.2%) [ig, S'/z/M^dm is a 
reasonable approximation for the final spin. One might 
expect that the conformal Kerr extrinsic curvature leads 
to the largest intrinsic spin, but as shown in Fig. [TJ the 
actual minimum in the ADM mass occurs at A = 0.7831, 
or roughly a mixture of 78% conformal Kerr and 22% 
BY extrinsic curvatures. As seen in Fig. [2J if we take 
A — 0.7831 and vary the mass parameter M^ in the ex- 
trinsic curvature (note that this is not the mass of the 
BH, rather it is the proportionality factor between S^nvc 
and the Kerr a parameter in the extrinsic curvature SKij) 
we find that the ADM mass is minimized when Mk — M . 
On the other hand, if we take A = 1, the ADM mass is 
minimized at slightly larger values of Mk. Finally, we 
need not have the spin of the wave (i.e. the spin param- 
eter in SKij) match the BY spin. In Fig. k3^ we plot the 
ADM mass for a BH of spin S = M^ and vary the "spin" 
of the curvature wave. The minimum in the ADM mass 
occurs when Swavc — S when A — 0.7831 and at slightly 
smaller values for A = 1 (note that A = 1 only correspond 
to pure conformal Kerr data when S'wavo = S) . 



The smallest ADM mass, and therefore the largest in- 
trinsic spins, are obtained when S'wavo = S, Mk = M, 
and A — 0.7831. The corresponding spin is Sn/M^^y^ — 
0.9324, which is slightly larger than the maximum spin 
of 0.932 obtained from purely conformal Kerr extrinsic 
curvature. 
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FIG. 1: The ADM mass of a single spinning BH with spin pa- 
rameter S — M^ versus weight A. Here the spin of the wave 
5Kij is fixed at S'wavo = M"^ . Note that A = 1, which corre- 
sponds to conformal Kerr extrinsic curvature, does not lead 
to the minimum mass, and hence maximum intrinsic spin. 
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FIG. 3: The ADM mass of a spinning BH with spin parameter 
S = A/^ as a function of spin of the curvature wave Swave for 
A = 0.7831 and A = 1. Here Mk = M. 



TABLE I: Initial spin parameters for Bowen-York data (BY), 
the new data at A™^", and conformal Kerr data. Note that 
^ — x/(l + \/l — X^)- In S'U cases, the spin SjM^ is identically 
1. 



ID 

B^^ 

A max 

cKerr 



eg = g/MJDM 
0.9282 
0.9324 
0.93207 



X 

0.9831 
0.9856 
0.9854 



SIhP 



Ho_ 



g = g/(2Mi 
0.8311 
0.8431 
0.8421 
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FIG. 2: The ADM mass of a single spinning BH with spin 
parameter S = M^ as a function of the Kerr mass parameter 
in the extrinsic curvature wave. Here the spin of the wave 
SKij is fixed at 5wavo ~ M'^. A minimum in the ADM mass 
is obtained for Mk ~ 1 when A = 0.7831. The minimum is at 
slightly larger values of Mk for pure conformal Kerr extrinsic 
curvature. 



B. Horizon mass dependence 

The difference between the ADM mass and the initial 
horizon mass of the black holes provides a measure of the 
energy initially lying outside the horizon. This energy 
can potentially escape to infinity or be absorbed by the 
black hole. We thus expect that the evolution of these 
spinning black holes will lead to a stationary black hole 
with 



S 



^S initial 






*^ Xfinal *^ Xinitial 



s 



Ml 

tin 



(20) 



For the A'^'^'^ = 0.7831 case, 0.9324 < xfinai < 0.9856. 
We also investigate the minimum of the horizon mass 
when A is varied but S and Mk are set M^ and Af, re- 
spectively. The results are shown in Fig. Ill The mini- 
mum horizon mass occurs when A"^'*'^ = 0.753, which is 
close to where the minimum in the ADM mass occurs. 
We summarize the initial parameters in Table [11 



V. EVOLUTIONS 

We evolve these single (distorted) BH data-sets us- 
ing the LazEv [27J implementation of the moving punc- 
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FIG. 4: The horizon mass of a spinning BH with spin param- 
eter S = M^ as a function of A, the mixing parameter for BY 
and conformal Kerr curvatures. 



ture approach [21 [3] with the conformal function W = 
exp(— 20) suggested by Ref. [55]. For the runs presented 
here, we use centered, eighth-order finite differencing in 
space f55^ and a fourth-order Runge Kutta time integra- 
tor. (Note that we do not upwind the advcction terms.) 

Our code uses the Cactus/EinsteinToolkit [21} [22] 
infrastructure. For these runs, we used the Carpet [M^ 
mesh refinement driver to provide fixed mesh refinement 
(since the BHs do not move across the grid); initially 
starting with 15 levels of refinement. When a larger AH 
forms, we reduce the numbers of levels from 15 to 11. Our 
base (coarsest) grid extended from — 400A/ to 400Af in all 
directions with a resolution oih = lOM/3. We used the it 
and z— reflection symmetries to reduce the computational 
domain by a factor of 4. The resolution on the finest grid 
was h = A-f/lOOO. We took a constant Courant factor of 
dt/h = 1/4 in all grids. We were able to achieve good 
spin preservation by choosing a grid structure where the 
finest grid is roughly twice as wide as the AH. Here we 
extended to finest grid to ±0.021Af . 

We found that using a lower-order dissipation operator, 
in this case fifth-order, gave sufficient accuracy, while, at 
the same time, reducing the number of required buffer 
points at the refinement boundary, as well as the cost in 
walltime of the dissipation step. At these boundaries, we 
used 16 buffer points. Standard ninth-order dissipation 
would require 20 buffer zones, while seventh-order would 
not require additional buffer zones, it proved to give a 
more noisy waveform than using the fifth-order dissipa- 
tion stencil. 

The initial data are axisymmetric, which we exploit in 
order to reduce the number of spectral collocation points 
to 140 X 140 X 4. We then use the full spectral expansion, 
rather than the much faster interpolation techniques, to 
transfer the initial data to the numerical grid. 

We obtain accurate, convergent waveforms and horizon 



TABLE II: Final spin parameters for Bowen-York data (BY), 
the new data at A™^", and conformal Kerr data. Note that 
Oh = Sh/MJj^, Orad = JADM/(Ai"ADM - 5Mrad)^. The no- 
tation of the form 1.23(4 ± 5) used below is shorthand for 
1.234 ±0.005. 



ID 


Oh 


Orad 


W-'SM,,^ 


BY 

A max 

cKerr 


0.9309(0 ± 2) 
0.935(2 ±1) 
0.9348(2 ± 1) 


0.9308(4 ±1) 
0.9351(3 ±1) 
0.9347(6 ± 2) 


1.47(6 ±5) 
1.49(1 ±6) 
1.49(0 ±7) 


ID 


Mh^ 


Sh 


Mi„ 


BY 

A max 

cKerr 


1.0365(1 ±3) 
1.0341(5 ±7) 
1.0343(1 ±2) 


1.0001(5 ±5) 
1.000(2 ±1) 
1.0000(8 ± 4) 


0.8563(6 ±3) 
0.850(9 ±1) 
0.85138(9 ±5) 



parameters by evolving this system in conjunction with a 
modified 1-f log lapse and a modified Gamma-driver shift 
condition [21 |30l [3T], and an initial lapse a{t = 0) = 
2/(1 + iPbl)- ^^^ lapse and shift are evolved with 

{dt-P'd,)a = ^2aK, (21a) 

dtP" = (3/4)f'' - 77/?" , (21b) 

where we use rj = 1 for all simulations presented below. 

The initial AH spherical and very small (see Fig. [5|, 
with coordinate radius r ~ 0.009. A larger AH forms 
after the BH absorbs the excess radiation on the grid, 
leading to an oblate spheroid AH with coordinate radius 
^ 0.3 — 0.5. We drop levels of refinement once this larger 
AH begins to equilibrate. 

We measure radiated energy in terms of the radia- 
tive Weyl scalar ■04, using the formulas provided in 
Refs. [321 133|. However, rather than using the full 1JJ4, 
we decompose it into £ and m modes and solve for the 
radiated linear momentum, dropping terms with ^ > 5. 
The formulas in Refs. [321 133] are valid at r = 00. We 
extract the radiated energy-momentum at finite (r = 
50Af , 60M, • • • , lOOM) radius and extrapolate to r = 00 
using both linear and quadratic extrapolations. We use 
the difference of these two extrapolations as a measure 
of the error. 

The initial 3-metric for these black holes is equivalent 
to the Kerr 3-metric (in quasi-isotropic coordinate) with 
a non-trivial axisymmetric distortion wave. This initial 
distortion contains energy that is partially radiated to 
infinity and partially absorbed by the black hole. From 
Table H] we see that in the extreme case 2.8% of the total 
mass lies outside the black hole (here S = 1) 



Ma 



DM 



M 



Ho 



1.0356 - 1.0073 
0.0283. 



(22) 



This difference, -Ewavo — A/adm — Mho, is the energy 
associated with the distortion wave. From Table [ill we 
see that 95% of the energy in the wave is eventually ab- 
sorbed by the black hole. Of the total ADM mass, only 



0.15% is actually radiated to infinity. This means that 
the black hole transitions from a nearly maximally spin- 
ning object to a submaximal Kerr black hole with spin 
0.935, as shown in the above Table |TTj This effect is also 
shown in Fig. [5J where, at around t — 10. 7M of the 
evolution, a new external AH appears suddenly with a 
larger area and much larger coordinate radius. Initially, 
only one trapped surface exists, with a very small co- 
ordinate radius of 0.0095M. As the evolution proceeds 
this horizon acquires an oblate form with the equato- 
rial radius nearly double the polar radius. The jagged 
features associated with the smaller trapped surface at 
later times are due to poor resolution of this very small 
object (which reduces in polar radius from ^ 0.009Af 
to < 0.002Af). At t = 10.7 M an outer AH forms with 
polar radius ~ 0.14M and equatorial radius 0.27M. At 
the same time a third, shrinking trapped surface also ap- 
pears that eventually meets the original smaller trapped 
surface, after t — 28.2M. These newly formed horizons 
are oblate. We can also observe notable quasinormal os- 
cillations soon after their formation. 

In order to understand the jump in the area in Fig. [5] 
when the initial distortion is absorbed by the rotating 
hole we note that 



Ah 
S 



Stt 



Stt 



i + v/T^x^ 



(23) 



where A^ is the area of the horizon. By keeping S fixed 
and allowing Mh to vary, we obtain 



dA 



H 



A 



H 



y/l^^ V Mh 



SM 



H 



(24) 



In our case this leads to a relative increase in the hori- 
zon area of 15% = 5.64 x 2.685% (using the numbers of 
Table [TT| in agreement with the 12.4 — 18.1% jump ob- 
served in the bottom plot of Fig. [5] (from 29.8 initially 
to 34 and then to 36.4) for the Ah- Another interesting 
aspect of these horizons is the jump in the deviation from 
nearly spherical to spheroidal. We define the deviation 
functions D to measure the non-spherical shape of the 
horizons. We define the coordinate dependent deviation 

^r V equatorial '"polar j/ v equatorial ' ?"polar/i wnerC T IS 

the coordinate radius, and coordinate independent devi- 
ation JJc = (Coquatorial~Cpolar)/(Cequatorial + Cpolar)j wherC 

c is the invariant measure of the circumference (see |34j 
for formulas relating D^ to the intrinsic spin a/M for 
Kerr). For a spherical horizon, Dc must be zero in all 
coordinate systems, while Dr can differ from zero, with 
the difference being a measure of how distorted the co- 
ordinate system is in the neighborhood of the AH. 

The initial horizon is highly spinning. Therefore, we 
would expect it to show the largest deviation from spher- 
ical. However, as seen in Fig. [6J this horizon is more 
spherical (i.e. \Dc\ is smaller) than the slower spinning 
equilibrium horizon (slower in the sense S/Mfj is smaller. 



S itself is constant). Also note that the distortion of the 
coordinates, as measured by the relative difference be- 
tween the Dc and D^. are much smaller for the equilibrium 
horizon despite the fact that this horizon is less spherical. 
We also observe a notable constancy of Dc of the smallest 
horizon, beginning from the initial time slice, through the 
formation of the outer apparent horizon, even until later 
times when the shrinking internal horizon meets the in- 
ner horizon at around t = 28. 2M. As expected the outer 
horizon is very distorted when it first forms and then set- 
tles to the Dc of a Kerr horizon with the corresponding 
spinning parameter a/M — 0.935. 




FIG. 5: (Top) The polar and equatorial radii of the three 
AHs found for the A = 0.7831 data. (Bottom) The area of 
the same three AHs. All in units of the ADM mass. The 
dotted vertical line at t = 10. 7M represents the time when 
an external trapped surface first forms. 

In Figs. [7] and [8] and Table |III[ we show the wave- 
form and radiated energy per £ mode for pure BY, K^^^ 
and cKerr extrinsic curvature data. Note how for all 
modes with I > 2, the cKerr data radiated the most, 
while pure BY data radiates the least. The effect, how- 
ever, is essentially in the higher £ modes. For the A™^^ 
data, the ADM mass is lower than for BY, and the sys- 
tem radiates more than BY. Both effects lead to larger 
final spins for the A™^'' data when compared to BY data. 



TABLE III: The energy radiated {SE/M), per £ mode, for the BY, A'"'"', and cKerr data. Note that the BY data radiates the 
most in the £ — 2 mode, but least in all other modes. The total radiated energy (SE/M) over all modes is (1.4799±0.0010) x 10"'^, 
(1.4978 ±0.0010) X 10"^ (1.4983 ± 0.0010) x 10"^ for BY, A™^", and cKerr, respectively. 



Mode 


BY 


A max 


cKerr 


(2,0) 


(1.416 ±0.001) X 10"'' 


(1.349 ±0.001) X 10-^ 


(1.318 ±0.001) X 10"^ 


(3,0) 


(6.01 ±0.01) X 10"^ 


(1.392 ±0.005) X 10"'' 


(1.672 ±0.007) X 10"* 


(4,0) 


(3.56 ±0.11) X 10-^ 


(8.85 ±2.8) X lO-'' 


(1.200 ±0.039) X 10-5 


(5,0) 


(3.12 ±0.30) X 10"^ 


(5.38 ±0.60) X 10"^ 


(7.25 ±0.86) X 10"'^ 


(6,0) 


(1.13 ±0.89) X 10-« 


(2.63 ±1.22) X 10"* 


(3.88 ±2.03) X 10-** 


(7,0) 


(0.3 ±1.0) X 10"^ 


(1.19 ±1.82) X 10"" 


(1.39 ±3.58) X 10"" 
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The A™^^ data, however, does not radiate more than the 
cKerr data. Thus there is a trade-off between the lower 
ADM mass of the A™'*'^ data and total mass radiated for 
the cKerr data. Shghtly smaller final spin could be ob- 
tained by optimizing the A parameter, but the effect will 
be quite small. Interestingly, the BY data radiates more 
in the £ — 2 mode than the other data. It is only after 
summing over all modes, that the other data sets end up 
radiated more. 



100.0 



FIG. 6: The deviation of the AH shape from a pure sphere 
D, for the A = 0.7831 data, as measured using the coordinate 
radii Dr = (rcquatoriai-rpoiar)/(f equatorial ±rpoiar) and the Co- 
ordinate independent measure using the equatorial and polar 

CirCUmterenCeS Uc ^ (^Cequatorial Cpolarj/tCequatorial ± Cpolarj- 

Note how the initial horizon starts out as a coordinate sphere, 
which rapidly becomes oblate. While the invariant circumfer- 
ences indicate a nearly constant deviation from spherical even 
at t=0. The dotted lines show the Kerr values for Dc for the 
initial intrinsic spin of xo = 0.9856 (upper) and the equilib- 
rium spin of Xoo ~ 0.9352 (lower). 
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VI. DISCUSSION 



In this paper we studied conformally flat initial data 
with an extrinsic curvature ansatz that interpolates be- 
tween the Bowen-York solution and the conformal Kerr 
solution. We found that, given our ansatz, the maximum 
intrinsic spin, normalized by the ADM mass, does not 
correspond to a pure conformal Kerr extrinsic curvature, 
but rather to a weighted sum of roughly 78% conformal 
Kerr and 22% Bowen-York extrinsic curvatures. While 
if we normalize the intrinsic spin with the horizon mass. 



the maximum occurs at roughly 75% conformal Kerr and 
25% Bowen-York extrinsic curvatures. This result can 
be interpreted as implying that, even if the conformal 
Kerr extrinsic curvature is the ideal choice for confor- 
mally Kerr initial data (3-metric), for conformally flat 
initial data, a mixture including Bowen-York compen- 
sates in part for the distortions produced by conformal 
flatness, leading to slightly larger intrinsic spins. 

Other numerical explorations are possible using 
Eq. (11) above for different choices of w. One could also 



try a variational approach, where an extremum in the 
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FIG. 7: The (^,0) modes of ^4 for the BY, A™'"', and cKerr 
data. The plot is arranged in two columns with even £ modes 
to the left and odd £ modes to the right. The plots are ar- 
ranged vertically with lower £ modes above higher £ modes. 
For the even £ modes, the odd part of ^4 is zero and not 
shown, while for odd £ modes, the even part of 7/14 is zero and 
not shown. Note how the different initial data give essentially 
the same £ = 2 signal, but differ quantitatively in the higher 
£ modes. 
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tem of five coupled partial differential equations. Our 
approach is rather to try to keep all the benefits of the 
simplicity of the Bowen-York solutions, i.e., an analytic 
confornial extrinsic curvature, Eqs. ( 15 )-( 19 ), with a very 



FIG. 8: The energy per {£, 0) mode for the BY, A" 
cKerr data. 
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simple form that is easy to incorporate in any numerical 
code that solves the Hamiltonian constraint for the con- 
formal factor. It is the intentions of the authors to deliver 
an open source code incorporating these new data to the 
Einstein Toolkit collaboration [25] . 

Ultimately we are interested in the evolution of these 
initial data and the ratio x = S/Mfj at late times, 
since this represents the actual spin of the black holes 
in a black-hole-binary simulation (that is, the individual 
black holes will equilibrate on much smaller timescales 
than the binary inspiral). Our numerical evolutions find 
that only an small fraction of the available energy initially 
outside the black hole radiates to infinity. Thus from the 
two initial data indicators, eg and x the former, which is 
normalized by the ADM mass is the most accurate. We 
thus conclude that our family of initial data is limited to 
represent black holes with intrinsic spins a"^'^^ < 0.9352. 

We also find an important jump in the apparent hori- 
zon location in numerical coordinates when radiation falls 
in and the hole settles from an extreme near-spherical 
horizon to a submaximal Kerr hole with a much larger 
coordinate size (an increase in coordinate radius of over a 
factor of 30). One needs to take into account these jumps 
both in coordinate size and physical mass when design- 
ing simulations and especially when measuring proper- 
ties such as the mass, spin, linear momentum 37;, and 
when tracking the motion of binary horizons in a full nu- 
merical simulation. More invariant measures of the hori- 
zon properties like its area, Figjsj and ratio of polar to 
equatorial circumferences, Fig. [Gjstill display an initial 
quasi-spherical (small distortion) AH, with a discontin- 
uous jump at t = 10. 7M to a larger, highly-distorted 
horizon that, after some oscillation, settles to the spher- 
ical deviations corresponding to that of a Kerr hole with 
a/AI = 0.935. The moving puncture approach also al- 
lows us to keep track of the internal horizon connecting 
the initial and final ones even after the formation of the 
latter at i = 10.7M. Although the details of this hori- 
zon transition are specific of the problem we studied, any 
other initial data set for multi black holes is expected to 
have a radiation content that is partially absorbed and 
partially radiated to infinity, and hence present similar 
qualitative early evolutions to those presented here. 



ADM mass can be found by varying the function w, but 
this is beyond the scope of the current paper. Note that 
Dain and Friedrich j35| have described the general form 
of the extrinsic curvature for asymptotically fiat initial 
data. 

Another example of conformally flat initial data with 
different extrinsic curvature is provided by the thin- 
sandwich approach [36j . but this requires solving a sys- 



VII. ACKNOWLEDGMENTS 

We thank S. Dain and O. Ortiz for valuable discussions 
on the subject of this paper. We gratefully acknowledge 
the NSF for financial support from Grants No. PHY- 
0722315, No. PHY-0653303, No. PHY-0714388, No. 
PHY-0722703, No. DMS-0820923, No. PHY-0929114, 
PHY-0969855, PHY-0903782, and No. CDT1028087; 
and NASA for financial support from NASA Grant No. 



07-ATFP07-0158. Computational resources were pro- 
vided by the Ranger cluster at TACC (Teragrid alloca- 



tion TG-PHY060027N) and by NewHorizons at RIT. 



[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005), gr- [19 
qc/0507014. [20 

[2] M. Campanelli, C. O. Lousto, P. Marronetti, and 

Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006), gr- [21 
qc/0511048. 
[3] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and [22 
J. van Meter, Phys. Rev. Lett. 96, 111102 (2006), gr- 
qc/0511103. [23 

[4] L. Gou, J. E. McCHntock, M. J. Reid, J. A. Orosz, J. F. 

Steiner, et al., Astrophys. J. 742, 85 (2011), 1106.3690. [24 

[5] C. O. Lousto and Y. Zlochower, Phys. Rev. Lett. 106, 

041101 (2011), 1009.0292. [25 

[6] H. Nakano, Y. Zlochower, C. O. Lousto, and M. Cam- 

paneUi, Phys. Rev. D84, 124006 (2011), 1108.4421. [26 

[7] U. Sperhake, V. Cardoso, C. D. Ott, E. Schnetter, and 

H. Witek, Phys. Rev. D84, 084038 (2011), 1105.5391. [27 

[8] G. Lovelace, M. Boyle, M. A. Scheel, and B. Szilagyi, 

Class. Quant. Grav. 29, 045003 (2012), 1110.2229. [28 

[9] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 

Rev. D74, 041501 (R) (2006), gr-qc/0604012. [29 

[10] C. O. Lousto and Y. Zlochower, Phys. Rev. Lett. 107, 

231102 (2011), 1108.2009. [30 

[11] C. O. Lousto, Y. Zlochower, M. Dotti, and M. Volonteri 

(2012), 1201.1923. 
[12] S. Komossa, Adv. Astron. 2012, 364973 (2012), [31 

1202.1977. 
[13] S. Brandt and B. Briigmann, Phys. Rev. Lett. 78, 3606 [32 

(1997), gr-qc/9703066. 
[14] A. Garat and R. H. Price, Phys. Rev. D61, 124011 [33 

(2000), gr-qc/0002013. 
[15] J. M. Bowen and J. W. York, Jr., Phys. Rev. D21, 2047 [34 

(1980). 
[16] G. B. Cook and J. York, James W., Phys. Rev. D41, [35 

1077 (1990). 
[17] R. Arnowitt, S. Deser, and C. W. Misner, in Gravtta- [36 
tion: An Introduction to Current Research, edited by 
L. Witten (John Wiley, New York, 1962), pp. 227-^265, [37 
gr-qc/0405109. 
[18] S. Dain, C. O. Lousto, and R. Takahashi, Phys. Rev. 
D65, 104038 (2002), gr-qc/0201062. 



page: 



page: 



G. B. Cook, Living Rev. Rel. 3, 5 (2000), gr-qc/0007085. 
M. Ansorg, B. Briigmann, and W. Tichy, Phys. Rev. 
D70, 064011 (2004), gr-qc/0404056. 
Cactus Computational Toolkit home 
http : //cactuscode . org. 
Einstein Toolkit home 

http : //einsteintoolkit . org. 

F. Loffler, J. Faber, E. Bentivegna, T. Bode, P. Diener, 
et al. (2011), 1111.3344. 

E. Schnetter, S. H. Hawley, and I. Hawkc, Class. Quan- 
tum Grav. 21, 1465 (2004), gr-qc/0310042. 
J. Thornburg, Class. Quant. Grav. 21, 743 (2004), gr- 
qc/0306056. 

O. Dreyer, B. Krishnan, D. Shoemaker, and E. Schnetter, 
Phys. Rev. D67, 024018 (2003), gr-qc/0206008. 
Y. Zlochower, J. G. Baker, M. Campanelli, and C. O. 
Lousto, Phys. Rev. D72, 024021 (2005), gr-qc/0505055. 
P. Marronetti, W. Tichy, B. Briigmann, J. Gonzalez, and 
U. Sperhake, Phys. Rev. D77, 064010 (2008), 0709.2160. 

C. O. Lousto and Y. Zlochower, Phys. Rev. D77, 024034 
(2008), 0711.1165. 

M. Alcubierre, B. Briigmann, P. Diener, M. Koppitz, 

D. PoUney, E. Seidel, and R. Takahashi, Phys. Rev. D67, 
084023 (2003), gr-qc/0206072. 

J. R. van Meter, J. G. Baker, M. Koppitz, and D.-L Choi, 

Phys. Rev. D73, 124011 (2006), gr-qc/0605030. 

M. Campanelli and C. O. Lousto, Phys. Rev. D59, 

124022 (1999), gr-qc/9811019. 

C. O. Lousto and Y. Zlochower, Phys. Rev. D76, 

041502(R) (2007), gr-qc/0703061. 

M. Alcubierre et al., Phys. Rev. D72, 044004 (2005), 

gr-qc/0411149. 

S. Dain and H. Friedrich, Commun. Math. Phys. 222, 

569 (2001), gr-qc/0102047. 

H. P. Pfeiffer and J. York, James W., Phys. Rev. Lett. 

95, 091101 (2005), gr-qc/0504142. 

B. Krishnan, C. O. Lousto, and Y. Zlochower, Phys. Rev. 

D76, 081501 (2007), 0707.0876. 



